Calculations of ^He+p Elastic Cross Sections Using Microscopic Optical Potential 
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An approach to calculate microscopic optical potential (OP) with the real part obtained by a 
folding procedure and with the imaginary part inherent in the high-energy approximation (HEA) 
is applied to study the ^He+p elastic scattering data at energies of tens of MeV/nucleon (MeV/N). 
The neutron and proton density distributions obtained in different models for *He are utilized in the 
calculations of the differential cross sections. The role of the spin-orbit potential is studied. Com- 
parison of the calculations with the available experimental data on the elastic scattering differential 
cross sections at beam energies of 15.7, 26.25, 32, 66 and 73 MeV/N is performed. The problem of 
the ambiguities of the depths of each component of the optical potential is considered by means of 
the imposed physical criterion related to the known behavior of the volume integrals as functions of 
the incident energy. It is shown also that the role of the surface absorption is rather important, in 
particular for the lowest incident energies (e.g., 15.7 and 26.25 MeV/nucleon). 
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I. INTRODUCTION 

The experiments with intensive secondary radioactive 
nuclear beams have made it possible to investigate the 
structure of light nuclei near the neutron and proton drip 
lines as well as the mechanism of scattering of the weakly 
bound nuclei. A special attention has been paid to the 
neutron-rich isotopes of helium (^'^He), lithium (^^Li), 
berilium (^'*Be) and others, in which several neutrons 
are situated in the far extended nuclear periphery and 
form a " halo" . A widely used way to study the struc- 
ture of exotic nuclei is to analyze their elastic scattering 
on protons or nuclear targets at different energies. Here 
we would like to mention, for example, the experiments 
on scattering of helium isotopes on protons at incident 
energies Einc less than 100 MeV/N, namely, for ^He at 
energy 25.2lE [1 i i, ||, 38.3 |], 41.6 UM, i and 
71 MeV/N [13,0, for ®He at energy 15.7 n^25.2 0, 
32 0, [nl, 66 [10, im and 73 MeV/N d and 

also at energy 700 MeV/N for He and Li isotopes (e.g. 

[il[il[il[i3,[il). 

The experimental data on differential and total reac- 
tion cross sections of processes with light exotic nuclei 
have been analyzed using a variety of phenomenologi- 
cal and microscopic methods (e.g. 

M El, HI m S ia iS, [33, M M, 13 ■ Among the 
latter methods we note, e.g. the microscopic analysis 
based on the coordinate-space ^-matrix folding method 
H, m. Hi m, m, S as weU as works where 

the real part of OP is microscopically calculated u sing 
the folding^proach (e.g. [11 ll, [11 [11, [11, [13, IM 
[i2,[il[iini|)- Usually the imaginary part of the OP's 
and the spin-orbit (SO) terms have been determined phc- 
nomenologically. Thus, the OP's have a number of fit- 
ting parameters. For example, OP's have been used to 





elaborate the elastic differential cross sections of ^He-|-p, 
^Re+^Re {Eiric=25 MeV/N) [H and ^He+p and ^He-Hp 
{Einc < 100 MeV/N) [23[by means of the M3Y-Paris 
effective NN interaction j42l. l43l |4^ . In the calcula- 
tions the proton and neutron densities of the helium iso- 
topes obtained by Tanihata et al. and also in the 
Cluster-Orbital Shell-Model Approximation (COSMA) 
[M [m, [10, [H were applied. It was shown [23] that 
the elastic scattering is sensitive to different density dis- 
tributions used in the folding approach. 

In our previous work 40] in order to exclude the us- 
age of the phenomenological imaginary part of OP we 
have performed calculations of ^He-|-p elastic differential 
cross sections by means of the microscopic OP with the 
imaginary part taken from the OP derived in [H, [4^ on 
the basis of the HEA [s^, [HI [11. This method (Glauber 
approach) in its optical limit |52l ] makes it possible to 
obtain an analytic expression of the scattering amplitude 
with the eikonal phase in the form of the so-called profile 
function. The latter is proportional to the integral of the 
one-particle density distributions of the colliding systems, 
and the integration is performed along a straight-line tra- 
jectory of motion. Generally, the integral contains also 
the form factor of the NN scattering amplitude and thus 
its form is akin to that of the standard folding poten- 
tial with the NN potential instead of the NN amplitude. 
The NN amplitude itself is known from the experimental 
data and therefore, the usage of a profile function offers 
certain advantages over approaches based on the folding 
potential. So, in nuclear physics, the HEA amplitude 
is ap plie d to energies larger than hundred MeV/N (see, 
e-g- m [H [13l- However, in the last two decades the 
HEA was generalized and applied to lower energies. The 
prescription to calculate the profile function consists in 
a replacement of the straight-line trajectory impact pa- 
rameter b by the distance of closest approach Tc in the 
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Coulomb field or by the respective distance rc„ in the 
presence of the nuclear field (real part of OP). Doing so 
a reasonable agreement with the experimental data on 
the proton- and nucleus- nucleus reaction cross sections 
has been obtained in the region of energies from 10 to 
1000 MeV/N (see, e.g. [11 li, [si, HI, [53, S [Hi, 113] ■ 
However, this approach becomes fairly rough when one 
calculates differential cross sections and also the total 
cross sections at comparably low energies. Besides, in 
the case of the microscopic OP given in a form of tables, 
this approach needs a numerical solution of the classical 
equation of motion to get the corresponding trajectory of 
motion. For this reason the method is not efficient for ap- 
plications because of quite complicated calculations. In 
this case the better way is to explore the equivalent HEA 
optical potential outlined in [H, \^ and then to solve the 
respective Schrodinger equation numerically using a stan- 
dard code which enables one to get the exact scattering 
amplitude and the total reaction cross sections including 
the interference terms as well. 

We used this approach in ^Ad\ to get the microscopic 
HEA imaginary part of the OP (ImOP) and added the 
real part of OP (RcOP) |4l|,[^. The ReOP includes the 
direct term and the exchange one which involves non- 
linearity effects. Also, the role of the spin-orbit inter- 
action has been considered. Additionally, the density 
dependence of the effective NN interaction, as well as 
the sensitivity of the results to the predictions of dif- 
ferent theoretical models for the density of ^He have 
been studied. It was shown that the more sophisticated 
Large-Scale Sheh Model (LSSM) [3 [H density of ^He 
is the most preferable one because it has led to a bet- 
ter agreement with the data. It was concluded in [i^ 
that the use of the microscopic folding ReOP {V^) and 
the HEA ImOP (W^) has led to agreement with the data 
on ^He-|-p elastic scattering cross sections for 41.6 and 71 
MeV/N. However, the data at lowest energy 25.2 MeV/N 
have been explained only on a qualitative level which is 
related to the limitations of using the HEA ImOP for 
energies around and less than 25 MeV/N. This has led 
to the necessity to reduce strongly the depth of HEA 
ImOP. It was shown in [i^ that the OP in the form 
Uopt = NrV^ + iNiW" with both and cal- 
culated microscopically and using only two free parame- 
ters Nji and Nj which renormalize the ReOP and ImOP 
depths can be reasonably applied to calculations of scat- 
tering cross sections at energies i?mc < 100 MeV/N, such 
as 41.6 and 71 MeV/N. 

In the present work we apply the developed approach 
to study the existing experimental data on ^He-|-p elastic 
scattering cross sections at incident energies less than 100 
MeV/N. Various model densities of ^He, such as those 
obtained within the approach of Tanihata et al. [4^ . 
LSSM [4l,[4|] and the Jastrow correlation method (JCM) 
[111, [13] are used to calculate the OP's. We study the role 
of the smn-orbit terms and in addition to our previous 
study [40|, we consider two more parameters N^'^ and 
Nf^ (when necessary) which renormalize the depths of 



the real and imaginary parts of the SO potential, respec- 
tively. In addition, the nuclear surface effects are also 
studied by introducing an additional surface term in OP. 
This is related to investigations of the lowest energy limit 
of the applicability of the HEA OP in ^He-|-p elastic scat- 
tering. Also we pay attention to the energy dependence 
of the parameters and Nj as well as to the respec- 
tive volume integrals. We note the necessity to analyze 
the differential cross sections estimating simultaneously 
the values of the total reaction cross section. This would 
give an additional test of the various ingredients of the 
approach. 

The theoretical scheme to calculate microscopically the 
real and imaginary parts of the OP, as well as the spin- 
orbit term is given in Section [ill The results of the cal- 
culations of OP's and elastic scattering differential cross 
sections, including those from some methodical ones, and 
their discussion are given in Sec. IIIII The summary of 
the work and conclusions of the results are presented in 
Sec. [Wj 

II. THEORETICAL SCHEME 

A. Direct and exchange parts of the real OP 
(ReOP) 

Here we give briefly the main expressions for the real 
part of the nucleon-nucleus OP that is assumed to be a 
result of a single folding of the effective NN potential and 
the nuclear densities. It involves the direct and exchange 
parts (for more details, see, e.g. [4l|, [i^, [i^ and also 

my- 

l/^(r) = F^(r) + V^^(r). (1) 

In Eq. HI) the direct part {V^) is composed of the 
isoscalar (IS) and isovector (IV) contributrions, corre- 
spondingly: 

VPsir) = / P2(v2)9(E)F{p2)v^,{s)dv^, (2) 



VFv{r)= J Sp2{r2)9{E)F{p2)vli[{s)dr2, (3) 
where s = r + r2 , 

P2(r2) = P2,p(r2,p) + P2,ri(r2,n), (4) 
Sp2{r2) = P2,p(r2,p) - P2,n(r2,n)- (5) 

In Eqs. dH) and ((5|) P2,p(r2,p) and p2,n(i"2,n) are the proton 
and neutron densities of the target nucleus. In Eqs. ([2|) 
and Q the energy dependence of the effective NN inter- 
action is taken in the usually used form: 

g{E) = 1 - 0.003^. (6) 
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Also, for the NN potentials v^q and w,^ we use the ex- 
pression from [4^ for the CDM3Y6-type of the effective 
interaction based on the solution of the equation for the 
f/-matrix, in which the Paris NN potential has been used. 
The density dependence of the effective interaction is 
taken in the following form: 



F{p) = C 



1 



-/9p(r 



7p(r) 



(7) 



where C=0.2658, a=3.8033, /3=1.4099 fm^, and 7=4.0 
fm'^. 

The isoscalar part of the exchange contribution to the 
ReOP has the form (see, e.g. [40|): 



ViTir) - 9{E) J p2(r2,r2-s)F 



P2 



(-1) 



(8) 



where the density matrix p2 (v2 , V2 
imated by the expression: 



P2 (r2,r2-s 
with 



) - P2{ 



s) is usually approx- 



(9) 
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ji{x) — —ji{x) — — (sin a; — a; cos a:) 



(10) 



and Vqq (s) is the isoscalar part of the exchange contri- 
bution to the effective NN interaction. The local mo- 
mentum fc(r) of the incident nucleon in the field of the 
Coulomb and nuclear potential (ReOP) is [i^: 



k'ir) = ( I? ) [E..r 



Vc{r) - V{r)] 



1 



A, 



(11) 



Substituting Eq. pT|) in Eq. ([8]) the iteration procedure 
was used to get the final result for the folding poten- 
tial. One can see that in this procedure the required 
microscopic potential V{r) (that has to be calculated ac- 
cording to Eq. ([T])) appears in the expression for fc^(r) 
[Eq. pip ] and, correspondingly, in the integrand of the 
integral in Eq. ([5]), i.e. in the expression for the exchange 
contribution to the OP. Thus, non-linearity effects occur 
as typical ingredients of the model and they have to be 
taken carefully into account. In our consideration, for 
the highest energy 73 MeV/N eight iterations and for 
the lowest one 15.7 MeV/N thirteen iterations were large 
enough in the calculations of the folding potentials. 

In Eq. ^ kp,2 is the average relative momentum of a 
nucleon in a nucleus |63|, '■ 



5 

3p 



1/2 



(12) 



where we choose for the kinetic energy density r(p) the 
expression from the extended Thomas-Fermi approxima- 
tion 165;, ,66]: 



r{p) 



2/3 



[p,{r)] 



5/3 



|Vp,(r)p ^ V^p,(r) 
36p,(r) 3 



(13) 



valid for each kind of particles q = n,p. It is shown in 
[40| how the isovector part of the exchange ReOP can be 
obtained. 



B. Density distributions of ^He 

In the calculations of the OP's we use the following 
point-nuclcon density distributions of *He: 

i) the Tanihata densities deduced in [i^l by means of 
comparison of the measured total reaction cross section of 
^'*He-|--^^C at 800A MeV with the respective expression 
from [67j derived in the framework of the optical limit of 
the Glauber theory: 



X 

Ppoint 



1 

— exp 



1 (^ - 2) 
63 3 



exp 



(14) 



Here X = Z, N and the parameter values of a and b can 
be determined from 



2 *2 I 1 

a ~ a I 1 



6^ ^b*^ [1 



AT 



(15) 



where a*=1.53 fm and fe*=2.06 fm; hence a=1.43 fm and 
6=1.93 fm for ^He. So, the proton distribution is defined 
by the first term only, while an excess of neutrons is de- 
scribed by the additional second term. The rms radii of 
the point-proton and point-neutron densities of ^He are 
equal to 1.76 fm and 2.69 fm, correspondingly; 

ii) the LSSM densities calculated in a complex Ahio 
shell model space [13, El] using the Woods-Saxon (WS) 
basis of single-particle wave functions with realistic ex- 
ponential asymptotic behavior; 

iii) the densities obtained in [6ll . [6^ with accounting 
for the NN central-type short-range Jastrow correlations. 



C. Optical potential within the high-energy 
approximation 

In Ref. [13] the so-called complex HEA optical poten- 
tial has been applied to explain the available data on the 
^Hc+p elastic differential cross sections and energies less 
than 100 MeV/N. The HEA OP was derived in [il on the 
basis of the eikonal phase inherent in the optical limit of 
the Glauber theory. Then, by means of this potential or 
taking only its imaginary part together with the folding 
real part of OP, the cross sections were calculated us- 
ing the code DWUCK4 [6^ for solving the Schrodinger 
equation. Thus, we don't apply the Glauber theory for 
calculating the scattering amplitude at relatively low en- 
ergies but utilize the equivalent HEA OP to solve nu- 
merically the respective wave equation. In this case, the 
use of the ordinary Glauber theory leads to insuperable 
problems in performing integration in the eikonal phase 
mentioned in the Introduction. Indeed, there one should 
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take into account the distortion of the integration path 
along classical trajectories in the field of the Coulomb and 
nuclear potentials (see, e.g. [H, [Hi, [H, [sS, [sl, [sl, [sO] . 
At the same time, to calculate the HEA OP one can 
use the definition of the eikonal phase as an integral of 
the nucleon-nucleus potential over the trajectory of the 
straight-line propagation, and have to compare it with 
the corresponding Glauber expression for the phase in 
the optical limit approximation. Doing so, the HEA OP 
can be obtained as a folding of form factors of the nuclear 
density and the NN amplitude fNN{q) [11,111]: 



{2-kY 



X / dqq jo{qr)p2{q)fNN{q)- 
Jo 



(16) 



In pB]) and cxnn are, respectively, the NN total 

scattering cross section and the ratio of the real to imag- 
inary parts of the forward NN scattering amplitude both 
averaged over the isospin of the nucleus. They both have 
been parametrized in [stI [69j as functions of energies in 
a wide range from 10 MeV to 1 GeV and also at energies 
lower than 10 MeV. The values of these quantities can 
also account for the in-medium effect by a factor from 
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FIG. 1: Total ((a) and (a')), point-proton (b) and point- 
neutron (c) densities of *He from the model of Tanihata 
COSMA [13, [H, LSSM m,!!!] and JCM calculations 



D. The spin-orbit term 



Following Refs. [68|, [7l|, [73| the expression for the spin- 
orbit contribution to the OP can be written in the form: 



VLsir) = 2X1 



Vo 



I dfRjr) 
r dr 



l rf/j(0 
r dr 



(1-s), (17) 



where A^=2 fm^ is the squared pion Compton wave- 
length, Vb and Wq are the real and imaginary parts of 
the microscopic OP at r=0, and /(r) is the form of the 
real [fR{r)\ and imaginary [//(r)] parts of the microscopic 
OP taken as WS forms f{r,RR,an) and f{r,Ri,ai). In 
our calculations the parameters (half-radius Rr{Ri) and 
diffuseness a/? (a/)) are obtained by fitting the WS poten- 
tial to the microscopically calculated real and imaginary 
contributions to the OP V{r) and W{r). 



III. RESULTS AND DISCUSSION 

In this Section we present the results of the calcula- 
tions of the microscopic OP's and the respective ^He-|-p 
elastic scattering differential cross sections at energies 
Einc < 100 MeV/N. In principle, the OP's do not con- 
tain free parameters, but they depend on the density 
distribution of the target nucleus. This allows one to 
test advanced theoretical methods that give predictions 
for the density distribution. In the case of *He we used 
the semi-empirical model of Tanihata , the large-scale 



shell model [44 , l45l | , as well as the results of the approach 
(6ll[6^ within the JCM. In Fig. [1] in logarithmic and lin- 
ear scales are shown the proton Pp{r), neutron Pn{r) and 
matter p(r) densities of ^He obtained in different mod- 
els. Also, for comparison, the known COSMA densities 
[Hi [Ml are presented. We note that among them only the 
LSSM densities have a realistic exponential asymptotics, 
whereas the others have a Gaussian one. The results for 
the JCM densities are given for the value of the corre- 
lation parameter (3—2.5 fm^^ in the Jastrow correlation 
factor 1 — e~'^ ^ , where r is the distance between neu- 
trons. It was shown in Refs. [6l|j [121 that the inclusion 
of this factor causes a slight increase of the density in 
the central part of the nucleus. Simultaneously, as can 
be seen in Fig. [21 this leads to a small decrease of the 
depth of the imaginary part of OP in comparison with 
the case of the Tanihata density. In the same Figure we 
show as examples the real and imaginary parts 
of the ^He-l-p OP's for energies 15.7, 32 and 73 MeV/N 
calculated using different densities. is calculated by 
a folding procedure and within the HEA (see Sec- 
tion |TT| . It is seen that the increase of the energy leads 
to reduced depths and slopes of ReOP and ImOP. 



We calculated the ^He-l-p elastic scattering differential 
cross sections utilizing the program DWUCK4 and 
using the microscopically obtained real and imagi- 
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FIG. 2: Microscopic real part (V-^) of OP ((a), (b) and (c)) 
and HEA imaginary part {W^) ((a'), (b') and (c')) calculated 
using the LSSM, Tanihata and JCM (/3 = 2.5 fm"^) densities 
of *He for energies i?=15.7 (solid lines), 32 (dashed lines) and 
73 MeV/N (dotted lines). 



nary contributions to the optical potential: 
Uopt{r) = NRV^{r) + iNiW"{r) + 



dr 



(18) 



r dr 



(l.s), 



where and are the depths of the SO optical po- 
tential obtained simultaneously with Rn^Rj) and aij(a/) 
from the approximation of the volume real and imaginary 
OP's by Woods-Saxon form. 

So, further in the present work we consider the set of 
the N coefficients as parameters to be found out from 
comparisons with the experimental data. We consider 
such a model as the appropriate physical basis, which 
constraints the fitting procedure by the established model 
forms of searching potentials. We would like to empha- 
size here that in our work we do not aim a perfect agree- 
ment with the empirical data. In this sense, the introduc- 
tion of the fitting parameters (A^'s) related to the depths 



of the different components of the OP's can be consid- 
ered as a way to introduce a quantitative measure of the 
deviations of the predictions of our approach from the 
reality (e.g. the differences of iV's from unity for given 
energies, as can be seen below). 

The discussion that follows is based on the fitting pro- 
cedure, where the additionally introduced strength pa- 
rameters iVfl, TV/, N^'^, Nj'~' are varied step by step. 
So, we start from the case Nr^Ni=1, N§^^=N^'^=0, 
then fit successively both coefficients Nji and Nj, and 
after that the values of N^'^ and Nf^. First, we give in 
Fig. [3] the results of our methodical calculations of the 
cross sections for different energies (15.7, 26, 32, 66 and 
73 MeV/N) using the densities of ®He from LSSM, Tani- 
hata and JCM approaches in the case when Nfi=Ni=l 
and N^'-'—Nj'^=0 (i.e. without spin-orbit interaction). 
It can be seen that the behavior of the cross sections for 
a given energy and interval of angles is weakly sensitive 
to the choice of the model for the density of ^He. In spite 
of this uncertainty we choose for the further applications 
the LSSM density since it has a realistic exponential be- 
havior in the peripheral region of the nucleus. 

The second methodical study is a test of the effect of 
Jastrow central short-range NN correlations on mecha- 
nism of the considered process of scattering. As known, 
the main parameter that governs the contribution of 
these correlations is f3, and we change it in wide lim- 
its from 2.5 fm^^ to 50 fm~^. It is seen in Fig. 0] that 
these changes result in an increase of the neutron density 
of about 2.5 times in the central part of ^He but this has 
no important effect on the calculated OP's and on the 
shape of the respective differential cross sections. There- 
fore, in the further calculations we do not account for the 
short-range correlation effects. 

Later, as a next step, we allow the "depth" of each of 
the parts of the OP in our semi-microscopic models 
to vary in order to find the optimal values of the param- 
eters Nji, TV/, N^'-' and Nf'-' by a fitting procedure to 
the available experimental data for the cross sections. In 
Fig. [5] we present the results of our calculations of *He-|-p 
elastic scattering cross sections for various energies and 
the LSSM density with the fitted values of the parameters 
Nji, Nj, N^'-' and Nf'-'. The values of these renormal- 
ization parameters are given in Table |T] together with the 
predicted total reaction cross sections. The results ob- 
tained using the values of the parameters from the first 
line of this Table for each energy are given by solid line in 
Fig. [5l while those from the second line for each energy 
are given by dashed line. 

As is known, however, the problem of the ambiguity 
of the values of the parameters arises when the fitting 
procedure is applied to a limited number of experimen- 
tal data. For instance, in the case of the LSSM den- 
sity, the values of iV/j=1.0 and 0.9, and correspondingly 
iV/=0.236 and 0.1 (with iV|O=0.107 and iV/o=0.040) 
lead to similar results in the case of 15.7 MeV/N. For 
E^32 MeV/N the results are similar when Nr^LO, 
iV/=0.374 and Nr=0A38, A^/=0.036; for E=66 MeV/N 
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FIG. 4: Point- neutron density of *He (a), and of OP 
((b) and (b')) for energy _E=32 MeV/N and *He+p elastic 
scattering cross sections (c) at energies i?=15.7, 32 and 73 
MeV/N calculated using JCM densities of *He for three values 
of the correlation parameter /3. 



TABLE L The renormalization parameters Nn, Nj, N'^'-^ and 
Nj''' obtained by fitting the experimental data in Fig. [5] in 
the case of LSSM density. The energies are in MeV/N and 
the total reaction cross sections an are in mb. 



the results are similar when iV/j=0.876, A^/=0.071 and 
Nr=0.854:, 7V/=0.086; for E=73 MeV/N they are sim- 
ilar when 7Vfl=0.875, iV/=0.020; Nr=0M9, iV/=0.010 
(with 7V|°=0.009 and A^^f o=0.002). Our calculations 
produce similar results when using the Tanihata density. 
We note that in some cases it has been enough to vary 
only the volume part of the OP, i.e. the values of the pa- 
rameters Nn and Nj without the spin-orbit parts of the 
OP. When all four parameters N are fitted the results for 
a given energy are similar, as already mentioned above. 
Thus, the problem to choose the most physical values of 
the parameters N arises. It is known that because the 
procedure of fitting belongs to the class of the ill-posed 
problems (see, e.g. Ref. [23|)) it is necessary to impose 
some physical constraints on the choice of the set of pa- 
rameters N. One of them is the total cross section of 
scattering and reaction. However, the corresponding val- 
ues are missing at the energy interval considered in our 
work. To our knowledge, the total reaction cross section 



E 


Nr 


N! 








15.7 


1.0 


0.236 








603.6 


15.7 


0.9 


0.1 


0.107 


0.040 


693 


26 


0.422 


0.104 


0.090 


0.010 


275.11 


26 


0.439 


0.144 


0.087 


0.023 


377.22 


32 


0.438 


0.036 


0.096 





71.9 


32 


1.0 


0.374 








419.5 


66 


0.876 


0.071 








55.7 


66 


0.854 


0.086 








65.9 


73 


0.875 


0.02 








1.48 


73 


0.869 


0.01 


0.010 


0.002 


1.22 



aR of ^He-fp process is known only at energy 670 MeV 
and it is about 200 mb [ll |. 

Another physical criterion that has to be imposed on 
the choice of the values of the parameters N is the be- 
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FIG. 5: The *He+p elastic scattering cross sections at differ- 
ent energies calculated using Uopt [Eq. (|18|l ] for various values 
of the renormalization parameters Nh, Ni, N'^''^ and Nf'~^ 
(presented in Table |T]| giving the best agreement with the 
data. The used density of *He is LSSM. Experimental data 
are taken for 15.7 [H , 26 0] , 32 [13, [H , 66 [13, [H and 73 

Mev/N [i3,[i3,[il. 



havior of the volume integrals [4l| 



Jv 



47r 
1 



drr^[NRV^{r)l 
^ ^ J drr^[NiW"{r)] 



(19) 
(20) 



as functions of the energy. 

It has been pointed out (see, e.g. Romanovsky et al. 
[t^ I and references therein) that the values of the volume 
integral Jy decrease with the increase of the energy in 
the interval < E < 100 MeV/N, while Jw is almost 
constant in the same interval. Imposing this behavior of 
Jy and Jw on our OP's (i.e. on their "depth" parame- 
ters Nr and Ni), we obtain by the fitting procedure the 
values of the parameters given in Table [TTl 

The results of the calculations of the cross sections are 
presented in Fig. [6] for the case of the LSSM density to- 
gether with the volume integrals Jy and Jw as functions 
of the energy. The results obtained using the values of the 



TABLE II: The parameters Nr, Ni, 7V|° and iVf °, the voL 
ume integrals Jv and Jw (in MeV.fm^) as functions of the 
energy E (in MeV/N) (selected in correspondence to the be- 
havior shown in Fig.[SIb,c)), and the total reaction cross sec- 
tions an (in mb) for the ^He-l-p scattering in the case of LSSM 
density (the results of the fit are shown in Fig. [6ja)). 



E N_ 



R 



Ni 



Jv Jw 



15.7 0.630 0.064 0.139 0.070 411.1 58.6 722.0 

15.7 0.630 0.052 0.166 0.057 411.1 47.6 701.2 

26 0.644 0.128 0.035 0.026 377.7 84.35 381.2 

32 0.648 0.120 0.062 0.022 358.3 69 302.7 

66 0.852 0.131 344.2 45 95.2 

73 0.869 0.090 0.004 330.0 29 60.9 

73 0.869 0.063 0.010 330.0 20.25 43.9 



parameters from the first line of Table HTl for the energies 
15.7 and 73 MeV/N are given by solid line in Fig. [H^a), 
while those from the second line for these energies are 
given by dashed line. In comparison to the data in Ta- 
ble |T] one can see that the total reaction cross sections 
decrease monotonically with the energy increased. Also, 
we reach the smooth change of the values of the volume 
integrals with the energy increase. Moreover, with the 
energy increase one sees the monotonic increase of the 
renormalization coefficients Nn of the volume real part of 
OP together with an " average" decreasing of Nj inherent 
in the imaginary part of OP. Almost a regular behavior 
is obtained for the spin-orbit correction coefficients N^*^ 
and Nj'-' . So, the coefficient is going to 1 in coinci- 
dence with a general conception of a folding procedure. 
But the obtained small values of Nj and problems with 
the fitting of our OP at low energies deserves a special 
attention. 

It is known that the fitting of the phenomenological 
OP's to the data of proton scattering on light nuclei leads 
to "shallow" imaginary parts of the OP's whose depths 
are sufficiently smaller than that of the real part of the 
OP. This has been observed in our previous works (40l . 
[tsI for the case of ^He-|-p elastic scattering. This is the 
case also in our present calculations. Another remark is 
connected with the difficulties in the description of the 
cross sections at low energies. In this case we cannot fit 
the data using only the volume form of OP. Instead, if 
one adds the contribution of a surface part of OP, then 
a better agreement with the data can be achieved. We 
would like to remind that such an admixture had been 
used in the earlier applications of the phenomenological 
OP (see, e.g. (zl). 

For this reason we consider also the contribution of the 
surface potential: 

dV^(r) 

Kptir) = Uoptir) ~ liaNs^^, (21) 

where the first term in the right-hand side is the expres- 
sion for the OP given by Eq. (in which the ImOP 
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FIG. 6: The *He+p elastic scattering cross sections (a) at 
different energies using LSSM density of *He and parameters 
from Table HI] Experimental data are taken for 15.7 26 
[2, 32 [13, O, 66 [13, [U and 73 MeV/N [l^, [0, [iTTlie 
obtained values of the volume integrals Jv (b) and Jw (c) 
(given by points) are shown as functions of the incident en- 
ergy, while the dashed lines give the trend of this dependence. 
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is taken in the form of V^{r)) and the second term is 
responsible for the surface effects. We would like to note 
that, in particular, for the lowest incident energy, the 
combination of the microscopically folded real and imag- 
inary parts in the form of is more appropriate. In 
Eq. (j2T]) a is the diffuseness parameter of V^{r) fitted 
by WS form. 

We present in Fig.JTJ^a) our results for both the volume 
and all components of the imaginary part of the potential 
Ugpf{r) and in Fig.JTJJb) for the cross section in the case of 
E=15.7 MeV/N obtained using the LSSM density of **He. 
The calculations are performed by fitting the strength 
parameters Nr, Ni, iV|°, 7V/° entering Eqs. HI]) and 
d^T]) and the depth parameter Ns of the surface term of 
the OP [Eq. In this case Nr=1.078, A^/=0.036, 

Ar|O=ArSO^0, A^s=0.207, a=0.686 fm, (7^=791.1 mb. 
It is seen from Fig. [7] that the inclusion of the surface 
contribution to the imaginary part of the OP improves 
the agreement with the experimental data, especially for 
small angles and in the region of the cross section min- 
imum. Obviously, for more successful description of the 
cross sections at low energies (15.7 and 26 MeV/N) our 
method has to be modified and improved by an inclusion 
of virtual excitations of inelastic and decay channels of 
the reactions. 

From the results presented in this Section one can see 



FIG. 7: (a) Volume (dashed line) and total (solid line) imag- 
inary parts of the OP U'optif) [Eq. (|2ip ] calculated using 
the LSSM density of **He for energy £=15.7 MeV/N; (b) 
The *He-|-p elastic scattering cross section at energy i5=15.7 
MeV/N using LSSM density of **He with fitted values of Nr, 
Ni, Ni°, Nf° and Ns given in the text. Experimental data 
are taken from [l2i] . 



that a notable renormalization of the imaginary parts of 
the microscopic OP and the necessity of its shape correc- 
tion at lower energies are needed for a reliable explana- 
tion of the data. In this connection one should remind 
that both the folding and HEA potentials have the same 
physical origin, namely, they are one-particle folding po- 
tentials, and thus they do not account for more compli- 
cated dynamical processes. We have already mentioned 
the role of the inelastic and breakup channels. In the 
last years many works have appeared where amplitudes 
of these processes were calculated within the distorted 
wave approximation and also by using the coupled chan- 
nel methods. The latter provide a way of estimating the 
elastic scattering cross sections, too (see, e.g. [T^ and 
references therein). On the other hand, if one considers 
the elastic channel itself, the general and formally estab- 
lished concept of the Feshbach theory 77] can give us 
a basis for the following qualitative physical suggestion. 
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Indeed, in this theory the elastic scattering potential is 
composed of two parts, the bare potential composed from 
one-particle matrix elements and the so-called dynamical 
polarization potential. Then, transforming this concept 
onto our model of OP, one can suppose that the bare OP 
is the microscopically calculated Vsp — + iW^ po- 
tential having the strength Njij = 1. And the rest part 
Vpoi = Vfit — Vsp, being the difference between the fitted 
OP and Vsp may be identified with a polarization poten- 
tial. In the framework of this outline of the scattering 
mechanism one can compare, for example, the imaginary 
part Wpoi with the imaginary part Wbreakup obtained by 
fitting the breakup cross sections with the respective ex- 
perimental data (see, e.g. [iHl) to make conclusions on 
the contributions of the breakup channel to the whole 
picture of scattering. In fact, it is seen in Fig. [7] a broad 
minimum of the Imt/^p((r) around 7'=2.6 fm which il- 
lustrates qualitatively the strong effect of the breakup 
channel on the elastic scattering cross section. 

IV. SUMMARY AND CONCLUSIONS 

The results of the present work can be summarized as 
follows: 

i) The optical potentials and cross sections of *He+p 
elastic scattering were calculated at the energies of 15.7, 
26.25, 32, 66 and 73 MeV/N and comparison with the 
available experimental data was performed. 

(a) The direct and exchange parts of the real OP {V^) 
were calculated microscopically using the folding proce- 
dure and density dependent M3Y (CDM3Y6-type) effec- 
tive interaction based on the Paris NN potential. 

(b) The imaginary part of the OP {W^) was calculated 
using the high-energy approximation. 

(c) Three different model densities of protons and neu- 
trons in ^He were used in the calculations: the Tanihata 
densities [13 , the LSSM densitie s [4l . l4H and the densi- 
ties obtained in an approach [U |62| with accounting for 
the central-type NN short-range Jastrow correlations. 

(d) The spin-orbit contribution to the OP was also 
included in the calculations. 

(e) The ^He-|-p elastic scattering differential cross sec- 
tions and total reaction cross sections were calculated 
using the program DWUCK4 [esj . 

ii) The density and energy dependence of the effective 
NN interaction were studied. It was shown that the be- 
havior of the cross sections for given energies and interval 
of angles is weakly sensitive to the choice of the model 
for the ^He density. The further calculations of the cross 
sections were performed using the LSSM density since it 
has a realistic exponential behavior in a peripheral region 



of the nucleus. 

iii) It was shown that the effects of the Jastrow cen- 
tral short-range NN correlations on the OP's and on the 
shape of differential cross sections are weak. 

iv) We note that the regularization of the OP's used 
in this work by the introduction of the fitting parame- 
ters (iV's) can serve as a quantitative test of our method, 
but not as a tool to obtain a best agreement with the 
experimental data. The problem of the ambiguity of the 
values of the parameters iVfl, TV/, Nf^ , Nf^ (that give 
the "depth" of each component of the OP) when the fit- 
ting procedure is applied to a limited number of experi- 
mental data is considered. It was shown that, generally, 
at energies i? > 20 MeV/N a good agreement with the 
experimental data for the differential cross sections can 
be achieved by varying mainly the volume part of the 
OP neglecting the SO contribution. A physical criterion 
imposed in our work on the choice of the values of the 
parameters N was the known behavior (e.g. [ill) of the 
volume integrals Jv and Jw as functions of the incident 
energy in the interval < Emc < 100 MeV/N. Another 
criterion is related to the values of the total cross section 
of scattering and reaction. However, the corresponding 
empirical data for these values are missing at the energy 
interval considered in our work. 

v) It was shown that the difficulties arising in the ex- 
planation of the ^He-|-p cross sections at lower energies 
(e.g. 15.7 and 26.25 MeV/N) lead to the necessity to 
account for the effects of the nuclear surface (and, cor- 
respondingly, of the diffuse region of the OP). For this 
reason we included in the cross section calculations the 
surface component of the OP and applied it to the case of 
E—l'a.7 MeV/N. In our opinion, the account of the latter 
can be considered as an imitation of the breakup channel 
effects. A more successful explanation of the cross section 
at low energies could be given by inclusion of polarization 
contributions due to virtual excitations of inelastic and 
decay channels of the reactions. 
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